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PURPOSE 


The  purpose  of  this  research  is  to  furnish  knowledge  that  can  be 
used  in  the  development  of  objective  methods  of  analyzing  and  fore¬ 
casting  tropical  weather  and  in  understanding  the  dynamics  of  tropical 
weather  phenomena.  It  is  planned  that  the  objective  techniques  will  be 
designed  for  electronic  computation  in  order  to  gain  speed  and  accuracy, 
and  to  reduce  personnel  requirements  in  an  operational  situation.  The 
investigations  are  divided  into  the  following  tasks: 

(1)  Analyze  the  three-dimensional  structure,  the  changes 
with  time,  the  movement,  and  the  surface  weather,  of 
selected  cases  of  representative  meteorological 
phenomena,  by  using  conventional  methods.  If  prac¬ 
ticable,  Tiros  data  will  be  included  in  the  analysis. 

(2)  Investigate  and,  insofar  as  possible,  develop  objective 
analysis  techniques  applicable  to  tropical  phenomena 
and  compare  the  objective  analyses  with  the  analyses 

of  Task  (1)  above. 

(3)  Utilizing  the  analyses  of  Tasks  (1)  and  (2)  above, 
carry  out  dynamical  studies  of  such  topics  as  the 
forces  predominant  in  various  phenomena,  the  con¬ 
servation  of  fields  of  vorticity,  divergence  and 
deformation,  and  the  importance  of  orography  and 
low-level  energy  inputs. 
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ABSTRACT 


During  the  third  quarter,  analyses  of  the  objectively  computed 
layer  quantities  (relative  humidity,  stream  function,  divergence,  etc.) 
for  the  period  5-8  May  1959  were  completed.  The  degree  of  persistence 
of  the  stream  function  was  investigated  and  found  to  be  considerably 
greater  than  that  of  heights  during  this  period.  The  stream  function 
fields  were  also  compared  with  subjective  wind-contour  charts  inde¬ 
pendently  prepared  by  Dr.  Portig  of  the  University  of  Texas,  and  showed 
excellent  agreement.  The  vertical  motion  and  humidity  fields  were 
compared  with  surface  reports  of  clouds  and  precipitation,  and  were 
found  to  be  in  reasonably  good  agreement.  The  limited  coverage  provided 
by  surface  reports  is  a  definite  handicap  to  such  comparisons,  and  in¬ 
dicates  that  further  evaluations  should  be  made  with  satellite  data. 
Investigation  of  the  influence  of  scale  (size  of  computational  triangles) 
on  computed  values  of  divergence  and  vorticity  revealed  that  magnitude 

decreased,  on  the  average,  by  about  one-half,  as  area  increased  from 
4  5  2 

2  X  10  to  4  X  io  km  .  Beyond  this  size,  average  magnitudes  remained 

constant  in  an  unrealistic  fashion.  Thus  it  appears  that  an  upper  limit 

5  2 

of  size  is  approximately  4  X  10  km  ,  for  computations  of  this  type. 
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I  PUBLICATIONS,  LECTURES,  REPORTS,  AND  CONFERENCES 


Publications,  Lectures,  and  Reports 

An  abstract  entitled  "An  Objective  Stream  Function  for  Tropical 
Analysis,"  has  been  submitted  to  the  Third  Technical  Conference  on 
Hurricanes  and  Tropical  Meteorology  to  be  held  at  Mexico  City  in  June 
1963.  This  paper  will  review  the  stream  function  technique  developed 
in  this  contractual  work  and  will  generally  follow  the  material  presented 
in  the  Quarterly  Reports. 

Conferences 

During  the  second  Western  National  Meeting  of  the  AGU  (27-29 
December),  Lt .  Col.  James  E.  Sadler  of  the  National  Science  Foundation 
(Indian  Ocean  Expedition)  visited  the  Institute.  The  applications  of 
satellites  and  numerical  techniques  to  tropical  meteorology  were  dis¬ 
cussed. 

Mr.  Marvin  Lowenthal  of  the  U.S.  Army  Electronics  Research  and 
Development  Laboratory  visited  the  Institute  on  28  January  1963.  During 
his  visit,  the  work  performed  under  the  contract  was  reviewed  and 
evaluated,  and  general  plans  were  made  for  the  remainder  of  the  con¬ 
tractual  period. 
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II  FACTUAL  DATA 


A.  Task  1 

The  processing,  mapping,  and  analysis  of  radiosonde  data  for  the 
Caribbean  during  the  period  5-8  May  1959  have  been  completed  in  the 
form  of  layer  averages  of  height,  temperature,  relative  humidity,  and 
wind  components.  Several  examples  of  these  data  have  already  been 
illustrated  in  previous  Quarterly  Reports.  The  relative  humidity  has 
been  particularly  revealing,  and  the  heights  and  winds  have  been  valuable 
for  objectively  computing  the  stream  function  discussed  in  Task  2.  A 
more  detailed  discussion  of  these  layer-averaged  data  and  the  general 
atmospheric  motions  and  structure  for  this  period  will  be  deferred  to 
the  Final  Report. 

B.  Task  2 

1.  Stream-Function  Routine 

In  the  Second  Quarterly  Report  the  development  of  an  objective 
method  for  computing  a  stream  function  (|)  was  discussed  in  some  detail. 
The  technique  used  was  to  successively  correct  an  initial  guess  (ob¬ 
served  height  data)  in  accordance  with  the  equation 


"  ~4g  ~ 1-  [(vo  +  V  (Xo  '  xi)  '  (Uo  +  Ui>  (yo  "  yi>]  (1) 


where  u  and  v  are  wind  components,  f  is  the  Coriolis  parameter,  g  is 
gravity,  the  subscript  o  denotes  the  station  of  interest  and  i  denotes 
a  neighboring  station.  An  insight  into  the  workings  of  the  procedure 
may  be  gained  by  applying  Eq .  (1)  to  a  hypothetical  group  of  stations 
arranged  in  a  rectangular  grid.  Station  o  is  at  the  center,  station  1 
lies  a  distance  d  to  the  east,  station  2  lies  similarly  to  the  north, 

3  to  the  west,  and  4  to  the  south.  Equation  (1)  is  written  for  relating 
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station  o  with  each  of  the  four  neighbors  (i  =  1,  2,  3,  4) .  For  example 
for  i  =  1  and  f  considered  constant 


K  =  K  +  h  (vo  +  V  <_ 


d) 


(2) 


The  four  relations  for  i  =  1,  2,  3,  4  are  then  summed,  giving 


i=l  U 


2  r(v,  -  v3)  (u2  -  u4 


2d 


2d 


(3) 


The  quantity  in  brackets  is  a  finite-difference  expression  for  vorticity 
(obtained  from  observed  wind  components)  and  ^  ^  -  4  $o  is  a  finite- 
difference  form  of  the  Laplacian  for  the  stream  function.  Viewed  in 
this  way,  the  stream-function  procedure  utilizes  vorticity  given  by 
observed  winds,  along  with  values  of  f  given  initially  by  height  data, 
in  order  to  compute  This  usage  may  be  contrasted  to  that  in  extra- 

tropical  meteorology  where  Eq.  (3)  is  used  to  obtain  vorticity  from 
values  of  height . 

The  stream-function  procedure  is  fast,  fairly  insensitive  to 
data  distribution,  operates  quite  well  when  data  are  missing,  and  could 
easily  incorporate  aircraft,  constant-level  balloon,  or  satellite  wind 
data.  Stream-function  fields  for  the  Caribbean  area  have  been  computed 
for  the  period  5-8  May  1959. 

One  of  the  advantages  of  a  stream  function  having  the  units 
of  meters  is  that  it  permits  comparison  with  subjective  contour  analysis. 
Or,  Wilfried  H.  Portig  and  associates,  of  the  University  of  Texas,  who 
have  considerable  experience  in  tropical  meteorology,  have  prepared 
contour  maps  for  the  same  period  (5-8  May  1959) .  Comparisons  between 
these  independently  produced  analyses  and  the  stream-function  fields 
for  this  period  are  presently  being  conducted.  It  is  our  understanding 
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that  Dr.  Portig's  technique  is  based  primarily  on  winds  and  their  time 
continuity  and  to  a  lesser  degree  on  heights  corrected  by  use  of  thick¬ 
ness  charts.  His  analyses  for  the  775-mb  and  250-mb  levels — in  close 
but  not  identical  correspondence  to  the  Layers  II  and  V  of  Figs.  1  and 
3 — are  illustrated  in  Figs.  2  and  4.  As  can  be  seen  by  comparing  Fig. 

1  with  Fig.  2  and  Fig.  3  with  Fig.  4,  there  is  excellent  agreement, 
while  both  are  significantly  in  disagreement  with  the  observed  (un¬ 
corrected)  heights  shown  in  Figs.  5  and  6.  Other  comparisons  of  the 
two  methods  of  analysis  have  also  shown  excellent  agreement.  Any  dis¬ 
crepancies  noted  could  be  attributed  to  the  differences  between  layer 
and  constant- level  winds  or  to  the  subjectivity  involved  in  both  analyses 
over  areas  of  sparse  data. 

One  of  the  qualities  sought  in  a  meteorological  quantity  is 
persistence — that  is,  that  it  be  valid  for  a  certain  period  of  time. 

In  order  to  investigate  the  persistence  of  the  stream  function  and  of 
reported  heights,  the  root-mean- square  changes  and  correlations  over 
various  time  intervals  have  been  computed.  The  correlation  coefficient 
is  particularly  revealing  because  it  eliminates  consideration  of  the 
variations  of  the  mean  fields.  The  results  for  Layers  II  and  V  derived 
from  the  eight  synoptic  times  under  investigation  are  tabulated  in  Table 
I.  It  is  quite  apparent  from  this  table  that  the  stream  function 
possessed  considerably  greater  persistence  (correlation  coefficients 
between  0.72  and  0.96)  than  did  the  height  (correlation  coefficients 
between  0.72  to  0.86).  This  difference  in  persistence  was  particularly 
noticeable  in  Layer  V. 

2 .  Triangle  Routine 

With  this  triple-station  computation  routine,  fields  of 
vorticity,  divergence,  vertical  motion,  deformation,  and  geostrophic 
departures  have  been  prepared  for  the  period  5-8  May  1959.  The  results 
appear  very  promising  for  portraying  the  intrinsic  properties  of  the 
atmosphere.  The  vorticity  in  particular  has  shown  excellent  agreement 
with  the  stream-function  fields.  The  divergence  and  associated  vertical 
motions  are  more  difficult  to  assess  because  they  are  more  susceptible 
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FIG.  1  STREAM-FUNCTION  ANALYSIS  FOR  LAYER  II  (850-700  mb) 
AT  1200  GCT,  8  MAY  1959 


FIG.  2  WIND-CONTOUR  ANALYSIS  (Portly)  FOR  775-mb  LEVEL 


Table  I 


COMPARISON  OF  HEIGHT  (Z)  VERSUS  STREAM  FUNCTION  (£)  FOR  5-8  MAY  1959 

(30  Caribbean  Stations) 


Layer 

Period 

(hours) 

Root-Mean-Square  Change 
(meters) 

Correlation 

Z 

? 

z 

n 

II 

(850-700  mb) 

12 

13 

10 

0.80 

0.92 

24 

16 

16 

0.81 

0.87 

36 

23 

24 

0.75 

0.80 

48 

30 

33 

0.72 

0.72 

V 

(350-200  mb) 

12 

32 

19 

0.73 

0.96 

24 

40 

27 

0.86 

0.96 

36 

51 

35 

0.79 

0.94 

48 

55 

43 

0.80 

0.92 

to  erroneous  and  unrepresentative  winds.  The  positive  vertical  motions 
of  the  lower  layers  have  shown  general  agreement  with  bad  weather  and 
low  cloud  build-up,  although  positive  vertical  motions  when  in  combination 
with  very  low  humidities  had  correspondingly  good  weather.  One  must 
realize  that  vertical  motion  is  simply  one  factor  in  the  production  of 
weather;  humidity,  stability,  advection,  among  other  factors,  also  need 
to  be  considered.  The  computed  vertical  motions  at  700  mb  on  6  May 
1959  (1200  GCT)  are  shown  in  Fig.  7,  along  with  the  humidity  patterns 
of  Layer  II  and  the  surface  weather  information.  This  chart,  although 
the  worst  of  the  four-day  series,  still  shows  reasonable  agreement  be¬ 
tween  the  three  patterns.  An  interesting  feature  of  Fig.  7  is  that 
considerable  land-sea  activity  was  occurring  along  the  coastal  regions 
of  Mexico  and  Central  America.  Unfortunately,  these  observations  fell 
outside  the  radiosonde  network.  It  is  believed  that  a  reasonably 
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FIG.  7  VERTICAL  MOTIONS  (700  mb),  RELATIVE  HUMIDITY  (Layer  II) 
AND  SURFACE  WEATHER 


i 

f 

i 

V 
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definitive  evaluation  of  kinematic  vertical  motions  can  only  be  made  by 
comparing  them  to  satellite  cloud  photos  and  to  vertical  motions  computed 
by  other  methods.  It  is  hoped  that  these  comparisons  can  be  carried  out 
in  future  case  studies. 

C.  Task  3 

An  important  aspect  of  the  triangle  computations  (vorticity,  diver¬ 
gence,  etc.)  has  been  the  effect  of  the  varying  sizes  necessitated  by 
the  scarcity  and  unevenness  of  reporting  stations.  The  question  of 
magnitude  of  divergence  and  vorticity  versus  scale  is  not  only  important 
in  the  evaluation  of  the  triangle  routine,  but  has  considerable  bearing 
upon  the  extension  of  numerical  forecasting  techniques  to  the  tropical 

regions.  The  triangles  constructed  for  the  computations  varied  in  area 
4  6  2 

from  about  10  to  10  km  ,  although  39  of  the  43  triangles  were  smaller 
5  2 

than  4  X  io  km  .  The  absolute  values  of  the  vorticity  and  divergence 

for  each  of  these  triangles  have  been  averaged  over  the  eight  synoptic 

times  under  study  and  over  the  six  different  layers.  The  resulting 

distributions  for  both  the  vorticity  and  the  divergence  showed  a  distinct 

decrease  with  increasing  area  (approximately  proportional  to  the  inverse 

5  2 

of  an  effective  radius)  for  triangle  areas  less  than  about  4  X  10  km  . 

For  the  larger  triangles,  computed  magnitudes  of  divergence  and  vorticity 
appeared  to  be  constant  and  not  dependent  on  triangle  size.  This  inde¬ 
pendence  of  size  is  not  realistic,  for  reasons  which  will  be  presented 
in  detail  in  the  final  report.  At  this  time  we  will  merely  mention  that 

the  kinematical  technique  becomes  ineffective  for  triangles  having  areas 
5  2 

much  greater  than  10  km  . 

Average  magnitudes  of  vorticity,  divergence  and  deformation  have  been 
computed  for  the  entire  period,  5-8  May  1959,  for  each  of  the  six  layers 

(Table  II).  The  total  magnitudes  were  1.7  X  io  5  sec  1  for  vorticity, 

-5  -1  -5  -1 

0.9  X  io  sec  for  divergence,  and  1.2  -  1.5  X  io  sec  for  defor¬ 
mation,  the  ratio  of  vorticity  to  divergence  being  1.9. 

During  the  quarter  a  preliminary  experiment  was  carried  out  in  the 
use  of  the  streamfunction  routine  for  computing  a  velocity  potential, 
for  the  purpose  of  describing  the  divergent  wind  components.  If  wind 
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vectors  are  rotated  90°  counterclockwise  by  the  operation  Ik  x  ^  (where 
\V  =  (Be  X  Vf )  +  ^x)  one  obtains 


|k  X  y  =  (Jk  X  7x)  -  Vijr 


(4) 


The  computer  program  used  previously  to  obtain  \|r  can  now  be  used  to 
obtain  the  potential  function  x  (with  x  constant  over  the  field  as  a 
first  guess) .  The  results  indicate  that  the  rotational  wind  components 
are  suppressed,  but  that  the  divergent  and  translational  components  re¬ 
main.  It  is  believed  that  the  translation  property  can  be  omitted  by 
further  experimentation. 
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Ill  PROGRAM  FOR  THE  NEXT  QUARTER 


The  evaluation  of  the  triangle  routine  and  the  investigation  of  the 
effects  of  scale  will  be  continued.  Further  comparisons  of  objective 
and  subjective  analysis  are  planned.  As  a  result  of  letter  and  phone 
communications  with  Dr.  Portig,  we  expect  that  he  will  visit  the 
Institute,  if  feasible,  for  thorough  discussions  of  the  two  types  of 
analyses . 
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IV  PERSONNEL 


Name 


Time  Devoted  to  Project 


R.  M.  Endlich,  Project  Leader 

R.  L.  Mancuso, 

Research  Meteorologist 

J.  Weaver,  Meteorological  Aide 


Approximately  100  hours 

Approximately  325  hours 
Approximately  190  hours 


This  Quarter 
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